Evidence of Magnetically Driven Structural Phase Transition in Parent Compounds 
i?FeAsO (R — La, Sm, Gd, Tb): study of low-temperature X-ray diffraction 
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We report measurements of structural phase transition of four parent compounds i?FeAsO (R = 
La, Sm, Gd, and Tb) by means of low-temperature X-ray diffraction (XRD). Magnetic transition 
temperatures associated with Fe ions (Tjvi) are also determined from the temperature dependence 
of resistivity. As R is changed from La, through Sm and Gd, to Tb, both the c-axis and a-axis lattice 
constants decrease significantly. Meanwhile both the structural phase transition temperature (Ts) 
and Tjvi decrease monotonously. It is also found that the temperature gap between Ts and Tjvi 
becomes smaller when the distance between FeAs layer becomes shorter. This result is consistent 
with magnetically driven structural phase transition and suggests that the dimensionality have an 
important effect on the AFM ordering. 

PACS numbers: 78.70.Ck; 74.70.Dd; 74.62.Bf; 74.25.Ha 



I. INTRODUCTION 



The recent discovery of superconductivity in layered 
pnictide-oxide quaternary compounds ROTmPn (R = 
lanthinides, Tm = Fe, Ni, Pn = P, As) has sparked 
enormous interest in this class of materials 1 J 2 -i 3 -i 4 ^'£' 7 -i£. 
Besides this 1111-type layered compounds, superconduc- 
tivity was subsequently discovered in other iron-based 
layered compounds with similar FeAs(Se) layers, i.e., 
122 systems^, 111 systems^, and 11 systems 1 ^. In 
all the FeAs-based parent compounds, there is a struc- 
tural phase transition in the temperature range 100- 
200 K, and a stripe- type antifcrromagnctic (AFM) or- 
dering associated with Fe ions accompanying the struc- 
tural transitio n 12 ' 13 : 14 . Various chemical doping ap- 
proaches can suppress the structural transition and AFM 
order, and high-T c superconductivity consequently ap- 
pears. Meanwhile, low-T c superconductivity has been 
observed in undoped FcP-based 1 - and NiAs(P)-based&i£ 
compounds with similar layered structure, but there is 
neither structural transition nor AFM ordering associ- 
ated with Fe(Ni) ions in these compounds. This result 
implies that there is a relationship between structural 
transition/ AFM ordering and high-T c superconductiv- 
ity. Theoretically, the origin of the AFM order is still 
controversial. There are two different physical pictures. 
One suggests that the AFM order is a spin density wave 
(SDW) which is driven by Fermi surface nesting be- 
tween the electron pockets at M point and hole pock- 
ets at T point based on band structural calculation a 12 ' 16 . 
The other suggests that the AFM order stems from the 
short range magnetic exchange coupling between local 
moment a 17 ' 18 : 19 ' 20 . However, regardless of the origin of 
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the AFM ordering, the theories suggest that the super- 
conductivity is tied to the magnetism in the FeAs-based 
materia la 16 ' 17 ' 20 ' 21 . The investigation on the structural 
properties and AFM ordering of the parent compounds 
can shed light on the mechanism of high-T c superconduc- 
tivity. 

The structural and magnetic transitions in the FeAs- 
based parent compounds are deeply connected. For the 
first discovered 1111 type systems, neutron scattering 
studies on the i?FeAsO (R = La, Ce, Nd, and Sm) sam- 
ples have found that the structural phase transition oc- 
curs first as temperature decreases, and then magnetic or- 
dering associated with Fe ions follow o 13 ' 22 , in contrast to 
the case in 122-typc systems where both structural tran- 
sition and AFM order occur at the same temperature^ 4 -. 
Recent report on isotope effect also shows positive iso- 
tope effect on both T c and AFM ordering temperature^ 3 .. 
Some theoretic studies proposed that the structural tran- 
sition is directly driven by the AFM orde r 17 ' 19 : 24 . In 
particular, a theory based on a Heisenberg-type local mo- 
ment exchange model suggests that the structural transi- 
tion can be driven by a nematic Ising magnetic order due 
to the presence of intrinsic magnetic frustration 1 ^. The 
theory^ 9 - predicts that the difference between the struc- 
tural and magnetic transition temperatures is controlled 
by the magnetic coupling between layers: the difference 
becomes larger when the coupling is weakened. Recent 
neutron experiment results in 4% Ni-dopcd BaFc2As2 
have supported this prediction 2 ^. 

In this paper, we report the investigation of structural 
phase transition detected by means of low-temperature 
X-ray diffraction (XRD) in the parent compounds 
LaFeAsO, SmFeAsO, GdFeAsO, and TbFeAsO. The 
AFM order temperatures associated with Fe ions (T/vi) 
and the AFM ordering temperatures (Tn 2 ) associated 
with the magnetic rare earths Sm, Gd and Tb are also 
obtained by measuring magnetic susceptibility and trans- 
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port properties. A systematic comparison of the struc- 
tural transition temperature (Tg) with Tjvi is made. As 
R is changed from La, through Sm and Gd, to Tb, 
both the c-axis and a-axis lattice constants decrease 
significantly. Meanwhile both T$ and Tjvi decrease 
monotonously. It is also found that the temperature gap 
(Tg — Jjyi) becomes smaller when the distance between 
FeAs layers becomes shorter. Therefore, our experimen- 
tal results provide concrete evidence supporting the the- 
ory proposed in RefJ^ and suggest the dimensionality 
may have important effect on the AFM ordering and on 
the superconductivity mechanism as well. 



II. EXPERIMENTAL 

The polycrystallinc i?FeAsO (R = La, Sm, Gd, Tb) 
samples were synthesized by solid state reaction in vac- 
uum using powders of RAs, R2O3 (for TbFeAsO, Tb40y 
was used instead), FeAs and Fe2As. RAs was presynthe- 
sized by reacting stoichiometric R pieces and As powders 
in evacuated quartz tubes at 1223 K for 24 hours. FeAs 
and Fe2As were prepared by reacting stoichiometric Fe 
powders and As powders at 1023 K for 20 hours. The 
powders of these intermediate materials were weighed 
according to the stoichiometric ratio of i?FeAsO respec- 
tively, and then thoroughly mixed in an agate mortar. 
The mixtures were pressed into pellets under a pressure 
of 4000 kg/cm 2 . All the processes were operated in a 
glove box filled with high-purity argon. The pellets were 
scaled in evacuated quartz tubes and heated uniformly 
at 1433-1453 K for 40 hours. 

The sample purity was first checked by measurements 
of powder X-ray diffraction (XRD) at room tempera- 
ture using a D/Max-rA diffractometer with Cu K Q radi- 
ation and a graphite monochromator. Low temperature 
X-ray diffraction (LTXRD) measurements for powder 
specimens were performed using the RINT 2500 system, 
Rigaku Co. An X-ray beam was generated by a rotating 
Cu anode. The specimens were cooled by a 4 He gas cir- 
culating cryo-cooler and can be cooled down to about 10 
K. The temperature stability is better than 0.1% during 
the LTXRD measurements. At several temperatures en- 
tire profiles of reflection peaks were measured with a step 
size of 0.01°and a step-counting time of 6 s and refined 
by the Rietveld method using the reported crystal struc- 
ture. For some reflection planes X-ray diffraction mea- 
surements with a step size of 0.005°and a step-counting 
time of 60 s were performed to accumulate more counts 
at certain temperatures. From the observed profile the 
d value of (220) peak, the integrated intensity (I.I.) and 
also the full-width-at-half-maximum (FWHM) were ob- 
tained. In these analysis the profile was fitted to a 
Pseudo-Voigt function. 

The electrical resistivity was measured by a standard 
four-terminal method. The samples for transport prop- 
erty measurements were cut into a thin bar. The tem- 
perature dependence of d.c. magnetic susceptibility was 
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FIG. 1: (Color online) (a) X-ray powder diffraction pattern 
at room temperature of i?FeAsO (R = La, Sm, Gd, Tb); (b) 
the variations of lattice constants a and c with the radius of 
R ions. The radius of R ions is taken from Ref>^. 



measured on a Quantum Design magnetic property mea- 
surement system (MPMS-5) under magnetic field of 1000 
Oe. 



III. RESULTS AND DISCUSSION 

Fig. 1(a) shows the room-temperature XRD patterns 
of i?FeAsO (R = La, Sm, Gd, Tb) samples and Fig. 1(b) 
shows the variations of lattice constants a and c with 
the radius of R ions. For all the four parent compounds, 
the XRD peaks can be well indexed based on a tetrag- 
onal cell with the space group of PA/nmm (No. 129), 
which indicates that the samples are in a uniform single 
phase without obvious trace of impurity phases. As R is 
changed from La to Sm, Gd, and then Tb, all the peaks 
shift to larger 2t9's significantly, implying a remarkable 
shrinkage of lattice in both a-axis and c-axis directions. 
This result is consistent with the fact that the radius of R 
ions decreases gradually as R goes from the light to heavy 
rare earth elements 2 ^. It can also be found from Fig. 1(b) 
that the c-axis shrinks slightly more quickly than the a 
axis does. 

Fig. 2 shows the temperature dependence of resistiv- 
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FIG. 2: (Color online) Temperature dependent resistivity of 
7?FeAsO (R = La, Sm, Gd, Tb). Inset: the derivative of 
resistivity(dp/dT') as a function of temperature. 



ity for the i?FeAsO samples. The resistivity starts to 
drop around 120-150 K. To show the drop position more 
precisely, we calculated the derivative dp/ dT versus T as 
shown in the inset of Fig. 2. The resistivity for LaFeAsO 
shows an upturn at low temperatures, but it remains 
metallic for _RFcAsO with R = Sm, Gd, and Tb. Similar 
resistivity has been observed in other 1111-type parent 
compounds with magnetic rare earth elements^. Such 
a difference in low temperature resistivity has not been 
well understood yet. We define the characterization tem- 
perature Tjvi as the peak position in the curves of dp/dT 
versus T. As shown in the inset, T/vi decreases signifi- 
cantly as R is changed from La to Sm, Gd, and Tb. Neu- 
tron studies have confirmed that the resistivity anomaly 
is caused by the structural phase transition and the fol- 
lowing formation of antifcrromagnctic SDW stated. Pre- 
vious reports have proposed that the peak position in 
dp/dT corresponds to the AFM ordering of Fe ions mo- 
ments rather than the structural phase transitio n 27 ' 28 . 
Indeed, the studies of neutron diffraction reported that 
the AFM ordering temperature of LaFeAsO is about 135 
K, about 20 K lower than the structural phase tran- 
sition temperature T$ of 158 Ki 3 -. The Tjvi value of 
LaOFeAs determined from the resistivity is 132 K, con- 
sistent with the AFM order temperature reported by the 
neutron diffraction. The AFM order temperature deter- 
mined from the measurements of Mossbauer spectroscopy 
and /^SR relaxation is also in agreement with 2jvi, the 
peak temperature in the dp/dT versus T curves^. Thus 
we can regard Tjvi as the characteristic temperature at 
which the magnetic moments of Fe ions become AFM or- 
dered. For R = Sm, Gd, and Tb, T N1 is 133 K, 128 K, 
and 122 K, respectively. 

In order to obtain more information about the mag- 
netism associated with Fe ions and R 3+ ions as R is mag- 
netic rare earth elements other than La, the magnetic sus- 
ceptibility of i?FcAsO (R = Sm, Gd, Tb) was measured 
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FIG. 3: (Color online) (a) Temperature dependence of mag- 
netic susceptibility measured under H of 1000 Oe for Sm- 
FeAsO, GdFeAsO and TbFeAsO. Inset: enlarged plot for low 
temperatures to show the AFM transitions of R i+ ions; (b) 
The plot of x" 1 versus T for GdFeAsO and TbFeAsO. The 
linear behavior for T > Tn2 means that x can be fitted by the 
Curie- Weiss law very well. Inset: Enlarged plot for SmFeAsO 
to show the kink around the AFM ordering temperature (TWi) 
of Fe ions. 



under the magnetic field of 1000 Oe, as shown in Fig. 3. 
From the temperature dependence of magnetic suscepti- 
bility shown in Fig. 3(a) and its inset, it can been found 
that the GdFeAsO and TbFeAsO have much larger mag- 
netic susceptibility compared to SmFeAsO because Gd 3+ 
and Tb 3+ ions have much larger magnetic moments. At 
low temperatures, clear phase transitions caused by the 
AFM ordering of the magnetic moments of R ions can be 
found at Tjv 2 . The AFM order temperature of R ions, 
i.e., T/V2 determined in our measurements, is 5.56 K, 4.11 
K, and 2.54 K for SmFeAsO, GdFeAsO, and TbFeAsO, 
respectively, consistent with previous reports ^ 29 ' 30 . For 
SmFeAsO, as shown in the inset of Fig. 3(b), a kink as- 
sociated with the AFM ordering of Fe ions can be ob- 
served around 134 K, consistent with the Tjvi value of 
133 K within the experimental error. For GdFeAsO and 
TbFeAsO, this kind of AFM order associated with Fe ions 
is buried in the large magnetic signals from the i? 3+ ions. 
The magnetic contributions from the i? 3+ ions obey the 
Curie- Weiss law very well. As shown in the Fig. 3(b), the 
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inverse of x increases strictly linearly with T as T > T^2 ■ 
For the SmFeAsO sample, the magnetic susceptibility 
does not exhibit the Cuire- Weiss behavior because the 
contribution from the Fe ions which is linearly depen- 
dent on temperature is comparable to the contribution 
from the Sm 3+ ions^i. By fitting the Curie- Weiss law, 
we obtained that the effective magnetic moments p e ff 
arc 8.0, and 9.7 fx B /f-U. for GdFeAsO and TbFeAsO re- 
spectively, which are consistent with the theoretical val- 
ues of magnetic moments of free Gd 3+ and Tb 3+ ions. If 
we subtract the Curie- Weiss term which should originate 
from the contributions of R 3+ ions, we can also find a 
slight drop in the subtracted term (x~Xcw) around T of 
120-140 K for GdFeAsO and TbFeAsO (not shown here), 
where xcw is the Cuire- Weiss fitting function. But it is 
very hard to distinguish whether such a drop in x — xcw 
occurs at the AFM order temperature of Fe ions or the 
structural phase transition temperature. Therefore, we 
will take T/vi as the transition temperature associated 
with the AFM order of Fe ions. 

As mentioned above, the neutron studies^ 3 - have re- 
vealed that a structural phase transition occurs just 
before the AFM ordering in LaFeAsO. It is generally 
believed that the structure of the parent compounds 
i?FcAsO transforms from tetragonal to orthorhombic 
when the temperature is lower than the structural tran- 
sition temperature Tg. Such a structure phase transition 
can be detected by the splitting of (220) peak in the low 
temperature x-ray diffraction. Fig. 4 shows the temper- 
ature dependence of the (220) peak d value for i?FeAsO 
samples. The inset shows the intensity of (220) peak be- 
fore and after the structural phase transition. It can be 
seen that the (220) peak splits into two peaks when the 
temperature is lower than the structural phase transition 
temperature. From the temperature dependence of (220) 
peak d-value, the structural phase transition temperature 
T$ can be easily determined as the splitting point of this 
peak. As R changes from La to Tb, the (220) peak d 
value decreases significantly, and T$ decreases as well. 

We summarize the variations of the structural phase 
transition temperature (Ts) and AFM ordering temper- 
ature (Tjvi) in Fig. 5. The structural and physical pa- 
rameters for these four parent compounds are also listed 
in Table I. As R is changed from La, through Sm and 
Gd, to Tb, not only T$ and Tjvi decrease significantly, 
but the temperature gap between Ts and Tjvi also be- 
comes smaller, i.e., the AFM transition occurs at the tem- 
perature closer to the structural phase transition tem- 
perature. Actually the structural phase transition and 
AFM ordering happen simultaneously in the more three- 
dimensional 122 parent compounds like BaFe2As2. 

The experimental results can be understood within the 
theory proposed in RefJ^ where an effective Heisenberg- 
type magnetic exchange model, the so called J\-J2-Jz 
model with Ji, J2 and J z being the in-plane nearest 
neighbor (NN), in-plane next nearest neighbor (NNN) and 
out-of-plane magnetic exchange couplings respectively, 
explains both the structural and magnetic transitions as- 
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FIG. 4: (Color online) The temperature dependence of the 
(220) peak d value for LaFeAsO (a), SmFeAsO (b), GdFeAsO 
(c) and TbFeAsO (d). The insets show the intensity of (220) 
peak before (blue lines) and after (red lines) structure phase 
transition. 
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FIG. 5: (Color online) Plot of structural phase transition 
temperature Ts, and magnetic transition temperature Tjvi 
associated to the AFM order of Fe ions, versus the c-axis for 
BFeAsO (R = La, Sm, Gd, Tb). Inset: plot of (T S -T N1 )/T N1 
versus J z /J2- The open circles denote the experimental data 
of (Ts — Tni)/Tmi, and the solid squares denote the theoret- 
ical values. See text for details. 



sociated with the FeAs layers. In this model, a collinear 
AFM ground state is obtained when Ji is less than 2J2, 
and a nematic Ising magnetic order transition that breaks 
lattice rotational symmetry takes place at a tempera- 
ture equal to or higher than the collinear AFM transi- 
tion temperature^. By including the lattice and spin 
coupling, the nematic order naturally produces an or- 
thorhombic lattice distortion. Therefore, the model cap- 
tures both structural and magnetic transitions and sug- 
gests the structural transition be driven magnetically. A 
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TABLE I: Structural and Physical parameters of _RFeAsO (R 
= La, Sm, Gd, and Tb) 

Sample r{R 3+ ) (A) a (A) c (A) T s (K) T N1 (K) T N2 (K) 



LaFeAsO 
SmFeAsO 
GdFeAsO 
TbFeAsO 



1.16 
1.08 
1.05 
1.04 



4.0349 8.7366 

3.9385 8.4941 

3.9151 8.4660 

3.8994 8.4029 



160 
144 
135 
126 



132 

133 5.56 

128 4.11 

122 2.54 



quantitative prediction of this model is that the difference 
between the structural and AFM transition temperature 
is determined by the ratio between J z and J2, i.e., J z j Ji- 
The AFM transition temperature, T N i oc J 2 / \og(J 2 / J z ). 
Using the calculated results in RefJ^ and assuming that 
J\ ~ Jii we find that the ratio J z / J2 for the four typical 
parent compounds (La, Sm, Gd, Tb)FeAsO are (4.10, 
15.2, 25.1, 40.0) xlO -4 . The out-of-plane magnetic ex- 
change coupling increases quickly as the lattice constant 
along c-axis deceases, and therefore the temperature dif- 
ference between Ts and Tjvi should decrease according 
to this model. The result is plotted in the inset of Fig. 5. 
The experimental data of (T5 — Tni)/Tni are in agree- 
ment with the theoretical calculations. By measuring the 
spin wave gap around the wavevector (0, it, 0) as shown 
in Rcfi2£ in future neutron scattering experiments, the 
value of J z can also be independently obtained. 



IV. CONCLUSION 

In summary, we have studied the structural phase tran- 
sition of four parent compounds i?FcAsO (i?=La, Sm, 
Gd, Tb) by measuring low-temperature X-ray diffrac- 
tions. As R is changed from La to Tb, T s and T an as well 
as the c axis decrease significantly. Furthermore, the tem- 
perature difference between Ts and Tjvi becomes smaller 
as the c axis becomes shorter. According to the theo- 
retical calculations proposed in RefJ^, the out-of-planc 
magnetic exchange coupling increases quickly with the 
decreasing c-axis lattice constant and therefore the tem- 
perature difference between Ts and Tmi is significantly 
influenced. The experimental data of (Ts — Tjvi)/Tjvi are 
in agreement with the theoretical calculations. This re- 
sult supports the theoretical proposal that the structural 
phase transition is driven by a nematic Ising magnetic or- 
der due to the presence of intrinsic magnetic frustration, 
and indicates that the dimensionality could have impor- 
tant effect on the AFM ordering, and the magnetism may 
play an important role in high-T c superconductivity of 
the iron pnictides. 
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